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Abstract We investigate the qubit geometric phase and its properties in de- 
pendence on the mechanism for decoherence of a qubit weakly coupled to 
its environment. We consider two sources of decoherence: dephasing coupling 
(without exchange of energy with environment) and dissipative coupling (with 
exchange of energy). Reduced dynamics of the qubit is studied in terms of the 
rigorous Davies Markovian quantum master equation, both at zero and non- 
zero temperature. For pure dephasing coupling, the geometric phase varies 
monotonically with respect to the polar angle (in the Bloch sphere representa- 
tion) parameterizing an initial state of the qubit. Moreover, it is antisymmetric 
about some points on the geometric phase-polar angle plane. This is in dis- 
tinct contrast to the case of dissipative coupling for which the variation of the 
geometric phase with respect to the polar angle typically is non-monotonic, 
displaying local extrema and is not antisymmetric. Sensitivity of the geometric 
phase to details of the decoherence source can make it a tool for testing the 
nature of the qubit-environment interaction. 

Keywords geometric phase • dephasing • dissipation • open system • Davies 
theory 

PACS 03.65.Vf ■ 03.65.Yz • 03.67.Lx 
1 Introduction 

One of the key obstructions of an effective implementation of quantum algo- 
rithms is related to the ubiquitous problem of decoherence in real quantum 
objects pp. Quantum decoherence is generic as it results from the imperfect 
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isolation of the quantum system from its environment. Decoherence can be 
diminished under very special conditions such as e.g. the presence of the de- 
coherence free subspaces [2] or via the application of tailored, external control 
schemes [3J. A promising novel direction in quantum information relates to 
so called holonomic or topological quantum computations [U[5] allowing for 
a substantial reduction of decoherence [51(7] • The essence of this method con- 
sists in encoding the information in the holonomy related to the geometric 
phase of the quantum evolution [8] . The geometric phase can be expressed as 
a path integral and via the Stokes theorem, can be converted into a surface 
integral. Therefore, it behaves like a geometric area. A quantity like an area is 
less dependent on the details of time evolution and therefore is less affected by 
changes of environmental conditions or an imperfect control, and hence, is typ- 
ically more robust. This is the key attribute that makes geometric phases at- 
tractive for the implementations of fault-tolerant quantum computation. Some 
suggestions have been presented to realize this objective, e.g. in NMR experi- 
ments [5] , ion traps [9] , neutral atoms in cavity QED [TU] , quantum dots [TT] 
or Josephson junction devices |12) . The performance of holonomic quantum 
gates under various conditions has been studied recently [T5] . 

The quantum evolution in the presence of decoherence is generically non- 
unitary. Therefore, the notion of geometric phase needs to be extended. There 
are several extensions of the geometric phase concept for systems which are 
either in a mixed state or/and undergo a non-unitary evolution. The first at- 
tempt towards this goal is given in [TJ] , being rather of mathematical charac- 
ter. The other are based on quantum trajectories |15j . quantum interferometry 
[IB] and the state purification (kinematic approach) [17iri8] . For non-unitary 
quantum evolution there is no commonly accepted scheme of defining the ge- 
ometric phase in open quantum systems [TJ5] . Here we use the approach based 
on state purification as proposed in Ref. [T5] ■ This so defined geometric phase 
has been extensively studied in various contexts [2"Dll2"T] . One of the appealing 
'advantages' of studying the phase defined in [18] is that it can be measured 
with a carefully prepared interferometric experiment [161118] . Our reasoning is 
thus guided by its potential for experimental implementation. 

There is no unique method of describing the time evolution of open quan- 
tum systems and there are several schemes to treat such systems which however 
typically give rise to non-equivalent dynamics [2"2"ll23| . One scheme consists in 
the derivation of a reduced system dynamics, via tracing over the degrees of 
freedom of the environment. Except some few exactly solvable models [H] 
|2~3] it is not clear how to relate the reduced dynamics to the microscopically 
first principles dynamics based upon the Hamiltonian structure of quantum 
dynamics [22,25 . The exactly solvable model of pure dephasing has been ap- 
plied for studying quantum channels [26] or in exploring the dynamics of quan- 
tum entanglement [2 7) . Despite its simplicity the highly non-trivial properties 
of geometric phase of qubits has also been discussed [§Ej. One of the most 
successful examples of constructing reduced dynamics is the Davies approxi- 
mation scheme |29j . Within this approach, starting with the general 'system- 
bath-interaction' Hamiltonian, one obtains, in a mathematically rigorous way, 
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a Markovian master equation form of a quantum system weakly coupled to 
the environment which preserves positivity and yields the correct equilibrium 
Gibbs state |30j . This approach has been applied to various problems of sta- 
tistical physics, quantum optics, solid state physics and quantum information, 
e. g. for studying entanglement dynamics in bipartite systems |31| . 

In this paper we apply the Davies master equation to study the geometric 
phase of the qubit coupled to the bosonic bath. Various families of a coupling 
and different coupling strength are shown to result in a qualitative and quan- 
titative modification of the geometric phase. This behavior could suggest a 
method to resolve the nature of the qubit-bath coupling: In particular, the 
dephasing coupling presents not only a mere theoretical construction but can 
be realized in experiments within tailored regimes |32j . In order to keep this 
study self-contained, we briefly review the notion of the geometric phase for 
a non-unitarily evolving qubit and then present the qubit master equation 
derived from the Davies theory. 



2 Geometric phase 

Generally, the time evolution of the qubit reduced density matrix p(t) is neither 
unitary nor Markovian [221125] , It is constructed as the mapping 

p{t) = A(t,t )p(t ) (1) 

obeying some properties depending on the specific circumstances and approx- 
imations such as e.g. the celebrated complete positivity condition [33]. In or- 
der to exploit the approach to the geometric phase based on state purifica- 
tion [TB] we have to present the reduced density matrix ([T]) in the spectral- 
decomposition form 

2 

p(t) = Y,M t )\ w *(t))(m(t)\, (2) 

8=1 

where Pi(t) and \wi(t)) are the eigenvalues and the eigenvectors of the matrix 
p(t), respectively. The geometric phase <l>(t) associated with such an evolution 
is defined as follows 1181: 



$(t) = arg 



E[k(°)^w] 1/2 (^(°)i^w) 



1=1 

■t 



x exp(- / (wi{s)\wi{s))ds) 



(i 



(3) 



where arg denotes the argument of the complex number, {wi\wj} is a scalar 
product and the dot indicates the derivative with respect to time s. For conve- 
nience, we assume the initial time being to = 0. For the sake of completeness 
we sketch here, following Ref. [15], the derivation of Eq.Q. The mixed state 
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defined by the density matrix ^ can be lifted to a pure state |<^(i)) in a larger 
Hilbert space, i.e., 

2 

!*(*)) = 53 v5^M*)>®l°i). ( 4 ) 

i=l 

where the vectors |aj) span the Hilbert space of an arbitrary ancilla. This 
is known as a purification of the density matrix p(t) in the sense that p(t) 
is a partial trace of the density matrix (<?(i)| over the ancilla Hilbert 

space. With the time evolution of the purified system one can associate the 
'Pancharatnam' relative phase 

a(t) = arg(!P(0)|ff(f)> (5) 

which contains both the gauge-dependent part (a dynamical phase) and a 
gauge-independent part. The central result of [T5] is to extract from Eq.©, by 
a proper choice of the 'parallel transport condition', the purification-independent 
part which can be termed a geometric phase because it is gauge invariant and 
reduces to the known results in the limit of an unitary evolution [341135] . The 
final result is then given by Eq. (J3J. 

As mentioned in the Introduction, this phase - contrary to other attempts 
of extending the notion of geometric phase for a non-unitary evolving quantum 
system - has a direct physical meaning as it can be measured via interfero- 
metric experiments [18] . i.e. one can construct the purification of the quantum 
system such that the relative phase reduces to the geometric phase (j3J 
after suitably defined 'compensating unitary' cutting of the dynamical part of 
the relative phase [18] , 



3 Weak coupling regime of qubit reduced dynamics 

The evolution operator A{t,to) defined by equation ([lj or its infinitesimal 
generator C defined by the equation 

j t p{t)=C{p{t)} (6) 

can be obtained in a few cases only; namely for stylized, exactly solvable models 
or in the limiting regimes such as the weak coupling limit or the singular 
coupling limit [33]. We consider a qubit coupled to a bosonic environment at 
temperature T. The Hamiltonian of such a system is chosen in the form 22 : 

H = Hq ®1 + T® H b + Hi ®Vb, (7) 

f-QO 

H B = / E(k)a*<(k)a(k)dk 7 (8) 
Jo 

P oo 

Vb= g{k)[a) (k) + a{k)]dk. (9) 
Jo 
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The operators a^(k) and a(k) denote the creation and annihilation boson op- 
erators, respectively. The qubit Hamiltonian and the interaction are assumed 
to take the form 

Hq = ~a z , Hi = hfi x a x + hfi z a z , (10) 

where oi are the Pauli operators, e is the qubit energy splitting and the dimen- 
sionless parameters [i x and \x z are coupling constants. Let us remark that if 
fj, x = the qubit energy operator Hq is an integral of motion, i.e. it commutes 
with the total Hamiltonian H leaving the expectation value of the correspond- 
ing energy observable unchanged. This situation defines the well known exactly 
solvable model of pure dephasing [21] • A non- vanishing fi x is then character- 
istic for exchange of energy and related dissipation processes. 

For an uncorrelated initial state p(0) (g) w{0) taken as a product of an 
arbitrary qubit density matrix p(0) and the equilibrium Gibbs state of the 
environment w((3) = exp(— /37Js)/Tr[exp(— (3Hb)} with /3 = l/ksT (ks is 
the Boltzmann constant), the Davies approximation for the Markovian kernel 
yields the following Markovian master equation pZMBTi] 

j t p{t)=C H {p{t)} + C R { P {t)}, (11) 

where the 'conservative' and 'dissipative' parts read as follows 

l 

C H {p{t)} = - l -[{H Q + Y, ^/w)4A0,P(*)L (12) 
fc,i=-i 



k,l=-l 

x ([A^(t),4j + [A H!i0 (t)4j), (13) 



where, see in Refs. [331131] . 

A kl =V k HiV u P±i = |±l)(±l|, 
O h = (A fc - \i)/h, A±i = ±e/2. (14) 

The states 1 1) and | — 1} denote the excited state and ground state of the qubit, 
respectively. The quantity c(uS) is the Fourier transform of the autocorrelation 
function of the bath operator Vb calculated in the Gibbs state w(j3) of the 
bath, namely, 



c(w)= / e- tut Tr[w(P)V B e ltHB/R V B e- ltHB/h }dt (15) 

J — OO 

and its Hilbcrt transform defines the function s(uj) in the following way 

*m = f r ^ 

2tt J^oo x-u 
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where P indicates the Cauchy principal value of the integral. 

In order to treat the complex qubit-environment interaction encoded in 
g{k) in © it is convenient to introduce the spectral density 

D(oj) = J dk\g(k)\ 2 5(oj(k) -w). (17) 

We further limit our consideration to the strictly Ohmic environment for which 
this spectral density is linear with respect to ui for small frequencies and ex- 
hibits an exponential cut-off frequency uj c , thereby exhibiting no non-physical 
ultraviolet divergences. Explicitly, this spectral density reads 

D{u) = —uexp(-ui/u; c ), (18) 

where the dimensionless parameter a characterizes the strength of the envi- 
ronmental influence on the qubit. Within this choice [55II5T] 

4»>-t(m^^ + ")«*-h/«J o» 

and s(w) is determined via the relation in (|16[) . 

In principle one can solve Eq. (Ilip using the Bloch vector formalism to 
obtain the coupled evolution equations for mean values (ak{t)),k = x,y,z to 
obtain the reduced density matrix as p(t) = (1/2)[1 + (a x (t))a x + {o y {t))<j y + 
{<J z {t)}cr z ]. This form allows to extract the spectral decomposition ([2} and the 
phase <&(t). Such an explicit form of the geometric phase result is, however, 
rather cumbersome without exhibiting much physical insight. We thus refrain 
from presenting such analytical details, but present here the full analysis of 
the geometric phase by numerical means. 



4 Analysis of geometric phase 

From Eqs ([I])-© it follows that in order to determine the geometric phase at 
arbitrary time t > 0, we must specify the initial state of the qubit. We consider 
the following class of initial states 

\6) = cos(0/2)|l) + sin(0/2)| - 1), (20) 

where is the polar angle in the Bloch sphere representation. The correspond- 
ing initial statistical operator p(0) takes the form 

P (o) = \e){e\. (2i) 

One of the eigenvalues of this operator is zero, say P2(0) = in Eq. ([3]), and 
it does not contribute to the geometric phase. This simplifies Eq. ([3]) in that 
only one term of the sum survives. The evolution of the freely evolving qubit, 
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Fig. 1 (Color online) Dependence of the geometric phase <P = $>[T) on the initial state of 
the qubit which is parameterized by the angle 8 on the Bloch sphere. The qubit is coupled 
to a purely dephasing Ohmic environment, i.e. fi x = in 1101 . The remaining parameters 
are: a = 1CT 2 , hui c /e = 10 2 and T = 0. 

with jjL x = n z = in (fTU|) . is cyclic with the time-period T = lufije and it 
acquires the geometric phase [35] 

<P Q = tt[1 + cos(0)], mod(27r), (22) 

which can serve as a reference for studying the influence of the environment. 
In the case of a coupling to an environment, the evolution of the qubit is not 
cyclic any longer. However, below we consider the phase = $(T) after the 
time T = 2Trh/e in order to study the role of coupling to the environment and 
for comparison with (|2"2"|) . 

The simplest situation occurs for pure dephasing; i.e. when fi x — so that 
the qubit energy does not change. The results presented in Fig. [T] show that 
the geometric phase plotted as a function of the initial state of the qubit (i.e. 
as a function of the parameter 9 in Eq. ({20)0 approaches zero (modulo 27r) 
with increasing coupling strength \i z . We observe that it varies drastically in 
the regime near 8 = n/2 and varies weakly outside this region. Moreover, the 
phase vanishes for — > (i.e. for the initially excited state |1) ) and 8 — > ir (i.e. 
for the ground state | — 1) ). This finding corroborates the results for the phase 
in the exactly solvable model of pure dephasing with arbitrary (not only weak) 
coupling [28]. For the presentation as in Fig. 1, we note that the function <P (8) 
is antisymmetric about the point {8,<P} — {ir/2,ir}, i.e. the relation 

$(- + 6) =2tt-$(^-8} for 8 e [0, tt/2] (23) 

holds. It can be interpreted as a rotation symmetry around the point {8, = 

W2.7T}. 




s 




Fig. 2 (Color online) Geometric phase •P = ^(T) vs. initial polar angle 9 for selected values 
of the dissipative qubit - Ohmic environment coupling strength /i x . The dephasing coupling 
strength is set at fi z = in HOP . and the remaining parameters are the same as in Fig. 1. 




Fig. 3 (Color online) Role of a qubit-environmcnt coupling on the geometric phase # = 
&(7~) vs. initial state preparation 8 for a qubit that is coupled to an Ohmic environment. 
The dissipative coupling strength is set at fi x = 0.3. The influence of dephasing is depicted 
for four different coupling strengths fi z . The remaining parameters are the same as in Fig. 
1. 



A most intriguing behavior on the role of the environment emerges when 
fj, x ^ 0; i.e. when the qubit-environment interaction is allowed to exchange en- 
ergy with the qubit system. The results depicted in Fig. [2] show the qualitative 
changes in the geometric phase properties for increasing dissipation coupling 
strength fi x as a function of the polar angle 6. Note that the geometric phase 
in Fig.[5]is plotted differently from Fig. [1] with <P/ir varying within the inter- 
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val (—1,1). We have decided to make this change in order to avoid confusing 
jump-like behavior of <P(6) in vicinity of the polar angle 9 = ir/2. E. g. the 
curve corresponding to the case \i x = 0.05 in Fig. [2] is very similar to the 
curve corresponding to the case \i z = 0.5 in Fig.[TJ However, presented in the 
'P/tt G (—1,1) interval it exhibits jump-like behavior which is an artefact of 
the way the plot is done. Let us recall again here that the phase <P is defined 
modulo 2ir. 

For small values of fi x , the geometric phase is close to that for the isolated 
qubit, cf. [i x = 0.05 in Fig. [2] when compared with Fig. [T]but with <P varying 
there <P/ir £ (0, 2). When /j, x increases the function <P(6) exhibits a local max- 
imum and minimum, see the case /i x = 0.3 in Fig. [2] For larger value of the 
coupling strength fi x (the case \i x = 0.4 ) the geometric phase is an increasing 
function of the polar angle 9 till to the value 9 = tt/2 reaching a local maxi- 
mum. Next, it decreases as 9 ~ > tt. In comparison to the dephasing coupling, 
in this case we can find at least three distinguishing features of the geomet- 
ric phase. Firstly, we note breaking of antisymmetry of < 1 > {9), being in distinct 
contrast to the case of pure dephasing (fi x — 0), cf. Fig. 1. Secondly, the depen- 
dence of the phase on the initial state parameterized by 9 is non-monotonic, 
exhibiting a local maximum and a minimum. Thirdly, the geometric phase <P 
vanishes for 9 — >• tt (i.e. for the ground state) but not necessary so for 9 — > 
(i.e. for the excited state). 

One can observe that for a fixed fj, X) the dephasing process controlled by 
IA Z does not change the qualitative properties of the geometric phase see 
Figs. 1 and 3. Pure dephasing affects only the off-diagonal elements of the 
density matrix, becoming closer to the maximally mixed qubit state. The 
geometric phase in a quantum evolution of such states vanishes. In a general 
energy-exchanging process the time dependence of the density matrix is more 
complex and the geometric phase <P seemingly quantifies this fact. Moreover, 
the stability of geometric phase with respect to decoherence is crucial for 
effectiveness of holonomic quantum computation |4llllj. It is evident that the 
stability of phase can be significantly improved via a proper choice of the initial 
state determined by 9 in Eq. (|2"0"j) . 

Thus far we considered zero temperature, T = 0. The effect of increasing 
temperature is depicted in Fig. 0] Firstly, we observe that if temperature in- 
creases the phase does not vanish for 9 — > tt while it tends to zero for 9 — > 0. 
Secondly, the main properties remain similar: In all presented cases a mini- 
mum and a maximum exist. However, the maximum diminishes with increasing 
temperature. 



5 Concluding remarks 

No realistic physical quantum system is in perfect isolation from its environ- 
ment. At best one can achieve a weak coupling between the system and the 
environment. In this weak coupling regime it is possible to extract the reduced 
dynamics of the open quantum system in a mathematically satisfactory and 
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Fig. 4 (Color online) The influence of varying temperature T on the geometric phase ^>(7~) 
vs. initial state preparation 9 is depicted for a dissipative qubit with fi x = 0.3 and zero 
dephasing, i.e. fi z = 0. The remaining parameters arc the same as in Fig. 1. Temperature is 
measured in units of e/ks- 



controlled way by using a Markovian reduced dynamics following the Davies 
scheme. In this work we have analyzed the geometric phase of a qubit in the 
presence of a weak coupling to a bosonic environment. We have investigated 
the relation between the geometric phase <!> and the mechanism for decoherence 
of the qubit for either the case of pure dephasing with fi x — or in presence 
of dissipative energy relaxation, i.e. fj, x ^ 0. The latter situation allows for a 
significant variation of the emerging geometric phase upon varying the cou- 
pling strength fi x . A variation of the pure dephasing coupling, i.e. fj, z ^ with 
fi x = 0, between qubit and environment barely affects the geometric phase. 
This feature is distinct from other set-ups, such as the emergence of quantum 
entanglement in open systems, where this dephasing-coupling mechanism can 
play a dominant or a similar role as an energy relaxation-coupling. 

Nowadays, the geometric phase plays a crucial role in a variety of physical 
problems and has observable consequences in a wide range of systems. Under 
various aspects, this concept occurs in geometry, astronomy, classical mechan- 
ics, and quantum theory. The impressive recent progress in nanotechnology 
and experimental techniques allows one to test the fundamentals of quantum 
dynamics and details of interactions modeled by Hamiltonians. The geometric 
phase is not a quantum mechanical observable, i.e. it is not represented by a 
Hermitian operator. However, it can be experimentally measured, cf. Ref. |36) . 
It can be used to encode information on systems. E.g. it has been proposed as 
an order parameter for quantum phase transitions [37] ■ The results obtained 
here suggest that one can also exploit the geometric phase as a quantifier char- 
acterizing a nature of the system-environment coupling. Indeed, three features 
of the geometric phase ^ = $(T) allow one to distinguish the character of 
qubit-environment coupling (i.e. pure dephasing vs. dissipation): (i) rotation 
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symmetry around some points on the 8 — <P plane or equivalently antisymmet- 
ric dependence of upon 8 about some points, (ii) non-monotonic behavior 
of $ with respect to 8 and (iii) the behavior of # for 8 — s> (i.e. for the qubit 
prepared in the excited state). We have verified that all these three features 
are manifest also at times t = nT(n — 2,3,4) for the measurable quantifier 
<P = <P(nT). This feature of the geometric phase <P thus presents an additional 
suitable tool in exploring characteristics of open system interactions at the 
quantum scale. 
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